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The study of randomness in low-dimensional quantum antiferromagnets is 
at the forefront of research in the field of strongly correlated electron sys- 
tems, yet there have been relatively few experimental model systems. Comple- 
mentary neutron scattering and numerical experiments demonstrate that the 
spin-diluted Heisenberg antiferromagnet La2Cui 2(Zn,Mg)204 is an excellent 
model material for square-lattice site percolation in the extreme quantum limit 
of spin one-half. Measurements of the ordered moment and spin correlations 
provide important quantitative information for tests of theories for this com- 
plex quantum-impurity problem. 
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The field of low-dimensional quantum magnetism has been of enormous interest to the 
condensed-matter physics community ever since the discovery that La2Cu04, the parent com- 
pound of the original high-temperature superconductor (La,Ba)2Cu04, is a model two-dimen- 
sional (2D) quantum (spin- 1/2) antiferromagnet. Because the superconductivity occurs in the 
vicinity of an antiferromagnetic phase in these materials, it appears likely that antiferromagnetic 
fluctuations are at least partially responsible for their rich physics. One of the new frontiers in 
condensed-matter physics lies in the field of quantum critical behavior, especially of "dirty" 
low-dimensional systems involving quantum impurities [|l]]. Although there has been much 
progress in the experimental investigation of quantum impurities in the simpler ID S = 1/2 
chain ^ and ladder compounds [|^], experiments with the 2D analog have been restricted to 
low impurity concentrations because of the lack of suitable samples [Q, ^, ^, 

We have investigated the properties of the spin- 1/2 square-lattice Heisenberg antiferromag- 
net (SLHAF) in the presence of a significant density z of quenched, spinless quantum impuri- 
ties, up to and through the percolation threshold. Specifically, the combined experimental and 
numerical results for the ordered moment Mst{z) and spin correlations ^(z, T) demonstrate that 
La2Cui_2(Zn,Mg)204 is well described by the Hamiltonian 

'H= J ^piPjSi-Sj, (1) 

{id) 

where the sum is over nearest-neighbor (NN) sites, J is the antiferromagnetic Cu-O-Cu su- 
perexchange, Sj is the S = 1/2 operator at site i, pi = 1 on magnetic sites, and = on 
non-magnetic sites. 

In the absence of quantum fluctuations, the NN square lattice undergoes a geometric transi- 
tion with site dilution z at the percolation threshold Zp 40.725% As indicated in Fig. 
1, below this concentration there is always one cluster of connected sites that spans the infinite 
lattice. Above Zp, the lattice consists entirely of finite-sized clusters. In studies of site-diluted 
5 = 5/2 Heisenberg and Ising antiferromagnets, long-range order was found to disappear only 
above Zp [[T0|]. However, in the extreme quantum limit oi S = l/2j)revious measurements of 
the Neel temperature T/v extrapolated to T/v = well below Zp [0, |^|^], and the possible exis- 
tence of a new quantum critical point at -25=1/2 < Zp has also been raised based on theoretical 
[ [TTi [T2| ] and numerical [ [T3| ] considerations. Recent Monte Carlo simulations of Eq. 1 suggest 



that the site-diluted 5 = 1/2 SLHAF remains ordered up to the percolation threshold [[Tj, ^5| ], 
but possibly with new, non-classical critical exponents [[T4|]. Our experimental data allow us 
to rule out the existence of a quantum critical point significantly below Zp, and we find that 
Mst{z) is reasonably well described, up to 2; ~ 30%, by a recent combined spin-wave theory 
and T-matrix approach [[T^ 1. 

In addition to information about the ground state properties, knowledge of the temperature 
and doping dependence of the 2D instantaneous spin-spin correlation length, ^{z, T), provides 
valuable information about the dynamics of the spin degrees of freedom. Neutron scattering 
studies of the spin correlations of La2Cu04 [ [T7| ] and related materials [ jlSl , [T^], as well as nu 



merical results [pO|], have served as important tests of theories for Eq. 1 in the absence of 



quantum impurities [|21], Remarkably, even though the quenched disorder leads to the loss 
of Lorentz invariance, a quantum non-linear sigma model approach developed for Eq. 1 in 
the absence of disorder is found to give an excellent effective description of ^(2;, T) up to at 
least z = 35%. Our experimental and numerical results can help guide the development of an 
underlying theory for the randomly diluted system. 

Several experimental difficulties have previously prevented quantitative experimental stud- 
ies of randomly diluted La2Cu04. First, the highest reported concentration of non-magnetic 
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ions is 25% [0, still well below Zp. Second, single crystal results have been limited even 
further, to 2; ^ 15% [^, 0]. Third, the excess oxygen typically found in as-grown samples intro- 
duces holes into the copper-oxygen sheets, which frustrate the antiferromagnetism and quickly 
destroy magnetic order [Z3Jl. Differing values for T]y{z) indicate that this problem has not been 
fully resolved [|, ||, §, |7|jr"^ 

By jointly substituting Zn and Mg on the Cu site, we were able to grow 
La2Cui_^(Zn,Mg)204+5 crystals by the traveling- solvent floating-zone method. The Zn con- 
tent was approximately 10% whereas the Mg content varied. Typical single grain sections were 
40 mm long and 4 mm in diameter. Both Zn^+ and Mg^+ are non-magnetic, effectively re- 
moving a magnetic site without introducing charge carriers. Zn^+ has a larger ionic radius than 
Cu^+, and Mg^+ has a smaller ionic radius than Cu^+. Compositional analysis was carried out 
by electron probe microanalysis on end sections of the crystals. The samples were carefully re- 
duced at T = 900 to 950 °C in Ar flow. Several crystals used for neutron scattering had a small 
Cu/Zn/Mg concentration gradient, typically 1 to 2% dilution, along the full length of the crys- 
tal. Mosaic widths were very good, 15' full width at half maximum (FWHM) or less. Small, 
very homogeneous sections a few cubic millimeters in size were cut from the larger crystals 
for magnetometry. Poly crystalline samples, with concentrations assumed to be equal to their 
nominal values, yielded similar magnetometry results. 

Elastic scattering at the (1,0,0) Bragg peak (orthorhombic notation), shown in Fig. 2A, 
was used to determine T/v and Mst(z). We find long-range Neel order even at z = 39%, where 
T/v is only 8 K. Figure 2B shows our result for T^^z). Below z ^ 20%, the data agree well 
with the linear behavior found in previous experimental ^, ^ and recent theoretical [ [T6| ] 
work. However, we find that T/v falls off more gradually at higher concentrations. Quantum 
fluctuations for S* = 1/2 apparently are not strong enough to noticeably shift the critical point: 
zs=i/2 = Zp, within the uncertainty of our experiment. 

In order to probe the ground state properties of La2Cui_2(Zn,Mg)^04 we have extracted 
Mst{z) from the low-temperature (1,0,0) intensity. To facilitate normalization between samples, 
we measured the intensity of a standard phonon for each sample to determine the illuminated 
volume; normalizing by sample mass gives the same result within the errors. The resulting 
normalized moment per Cu atom is plotted in Fig. 2C. Our data are consistent with previous 
nuclear quadrupole resonance(NQR) measurements below z = 15% [p. Numerical results 
for Eq. 1 deviate somewhat in their magnitude, but also suggest that -25=1/2 = Zp [ [T4l , p3| ]. 
Recent theory [|16p, which uses a spin-wave theory and T-matrix approach and is expected to 
be valid at low and intermediate concentrations, gives a good description of our data up to 
z ~ 30%. The classical (S 00) behavior is significantly different, because the reduction of 
the moment is attributed only to Cu sites being disconnected from the infinite cluster. It should 
be noted that already MstiO) is reduced from the classical value Mst,diO) = 1/2 by about 40% 
because of quantum fluctuations p^ , P4| ]. These fluctuations evidently increase with dilution, 
further weakening the ordered moment, but not sufficiently to disorder the system below Zp. 
The experimental data are well described by a power law with 2:5=1/2 = Zp and an effective 
exponent ^^.ff = 0.45(3). In their finite-size scaling analysis, Kato et al. [ [T4l ] extracted spin- 
dependent critical exponents, with P = 0.46 for S = 1/2, which is substantially different from 
the classical value Pci = 5/36 [^. However, it has been argued that these are not the true critical 
exponents, which should equal the classical values [[T5|]. 

We also performed a systematic study of the instantaneous 2D spin-spin correlation length 
^{z,T) in the paramagnetic phase of La2Cui_^(Zn,Mg)^Q4. The equal-time structure factor 
was measured in two-axis, energy-integrating mode ^1%. Figure 3A shows representative data 
for z = 19% at T = 200 K. The scattering broadens as ^ decreases, both with increasing 
temperature and increasing dilution. Correlation lengths were obtained from fits to a Lorentzian, 
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~ 1 /(I+O'Id^^)' where q2D is the 2D momentum transfer component in the Cu02 sheets relative 
to the zone center at i7 = 1, convoluted with the instrumental resolution. The extracted lengths 
are plotted versus J/T in Fig. 3B. The data shown were cut off above Tn by 1 SD 4 K), 
as obtained from fits of the order parameter assuming a Gaussian distribution in T/v (Fig. 2A). 
Temperature is scaled by J = 135 meV, the antiferromagnetic superexchange energy of the pure 
system [jT?]]. The data reveal that doping significantly decreases the rate at which correlations 
grow as tne system is cooled. Specifically, at high concentrations E,{z,T) crosses over from 
exponential to power-law behavior. We note that the z = 40(2)% and 43(2)% samples do 
not exhibit Neel order at 1.4 K (J/T ^ 1100), and the spin correlations appear to approach a 
constant zero-temperature value, as expected for z > Zp. 

In order to test the degree to which the experimental system is described by Eq. 1, we have 
performed quantum Monte Carlo (QMC) simulations to calculate ^{z, T). We used the loop- 



cluster algorithm, which is suited to study the randomly diluted Heisenberg model [12^]. Lattice 
sizes 10 to 20 times larger than the correlation length were used in order to avoid finite-size 
effects. We were able to perform our calculations on very large lattices of up to 1700 x 1700 
sites and to temperatures as low as T = J/100. Previous numerical work involved system sizes 
as large as 20 x 20 [|l3m] and reached T = J/2 [|2|] and T = J/20 [jTS]]. Between 5 and 
200 configurations were averaged for each temperature and concentration, and we chose 10^ to 
10^ equilibrations and lO'' to 10^ measurements per configuration. The QMC results (Fig. 3B) 
extend to higher temperatures and thus complement the experiment, covering a combined three 
orders of magnitude in temperature. We emphasize that this comparison contains no adjustable 
parameters, because J is known rather well for the experimental system. We find excellent 
quantitative agreement up to the percolation threshold. 

The ground state of the pure SLHAF is ordered, but quantum fluctuations renormalize the 

spin-wave velocity, c = 2\/2SZcJa, and spin-stiffness, ps = S'^ZpJ, from their classical values 
(we use units in which g^B = kB = h = 1). For S = 1/2 and z = 0, the quantum renormaliza- 
tion factors Zc and Z^ are known from various theoretical and numerical studies. Using these 
quantities, ^(z = 0, T) is quantitatively given by 



^ e c/a 
a 8 27cps 




(2) 



Even though Eq. 2 is strictly valid only at asymptotically low temperatures p^ , pj , E^, it 
agrees remarkably well with experiment [[T7|, [T3, [1^] and numerics [^] for 5* = 1/2 in the 
range 2 < ^/a < 200 shown in Fig. SB. The (derivation of Eq. 2 involves a mapping of the 
discrete Heisenberg Hamiltonian Eq. 1 to a quantum non-linear sigma model (QNLaM), and it 
is based on the assumptions of an ordered ground state and of translational invariance. Random 
dilution breaks translational invariance of the SLHAF and leads to defect rods in the time- 
like direction of the effective QNLcrM. The latter implies the loss of Lorentz invariance, and a 
QNLcrM description may no longer remain valid [|2T|, [T^. Nevertheless, it is valuable to test the 
extent to which this description may hold. Because there exist no accurate predictions for c{z) 
and Psiz), such a test requires us to treat these zero-temperature quantities as fit parameters. A 
modified form of Eq. 2, 

^ _ e c/a e^'^P'l^ 
a~ 8 27ipsl + {47ips/T)-'^ 

with ut = 1, has been suggested for disorder-free systems approaching a quantum critical point 
[^, where ps = and ^ ~ 1/T [pT]]. Even though zs=i/2 = Zp may not be a quantum 
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critical point, random dilution reduces the spin stiffness and may be viewed as bringing the 
system closer to such a point in an extended parameter space. In a classical picture, z = Zp is 
a (geometric and thermal) multi-critical point, and = as well as power-law behavior are 
expected. The thermal correlation length exponent ut = 0.90(5) was previously determined for 
the 5 = 5/2 system Rb2(Mn,Mg)F4 [jT^]. Furthermore, for z > Zp it was found that the data 
were well described by the form 

' ' (4) 



where is a zero-temperature length, possibly reduced from the geometric length of the per- 
colation problem by small quantum fluctuations, and ~ T^"^. 

We find that Eq. 2 gives a very good description of our data, especially at lower concentra- 
tions, and for ^/a > 8 to 10 at 2; = 31% and 35%. The modified form Eq. 3, however, even 
captures the high-temperature power-law behavior at the higher concentrations. Results of fits 
to Eq. 3 of our QMC data for z < Zp are shown in Fig. 3B and summarized in Table I. Fits 
to Eq. 2 result in large uncertainties in the value of the spin- wave velocity, but Ps(z) and c{z) 
extracted using these two forms agree within the errors. For the pure system, a fit of numerical 
data [|0|] below ^/a = 200 yields 27rp,(0) = 1.18(1) J and c(0) = 1.33(3) Ja, about 4% higher 
and 20% lower, respectively, than the most accurate estimate [^. A recent combined QNLaM, 
and percolation theory approach [ [T2| ] resulted in 



Ps{Oj 



m 



(5) 



where g{0) = 0.685 is the coupling constant corresponding to the S" = 1/2 SLHAF at 2 = 
and A{z) and Poo (-2) are the bond dilution factor and the probability of finding a spin in 
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le infinite cluster. The latter two quantities are well described up to 2 ^ 37% and z ^ 30%, 
respectively, by A(z)^ 1 ~ nz + nz"^ /2 and Poo{z) ~ 1 — z [|^, Equation 5 incorrectly 
predicts Ps(z'^ 30%) = 0, and hence a quantum critical point well below the percolation 
threshold [O]. To our surprise, we find that substituting 1 + z fox 1/Pqo{z) quantitatively 
describes p7(^)/p(0) even at 2; = 35%, as shown in Table 1. This substitution is correct at low 
concentrations and prevents the second term in Eq. 5 from going to zero below Zp. We note 
that Eq. 5 is a one-loop renormalization-group result, and higher-order terms can be expected 
to improve agreement with our observations. The expression for the spin-wave velocity that 
corresponds to Eq. 5 [[T^, c(2)/c(0) = A{z)(l + z/2), decreases monotonically with dilution, 
and hence does not describe the behavior found from our fits (see Table 1). 

For z > Zp, the joint experimental (z = 40(2)% and z = 43(2)%) and numerical {z = 41% 







and z = 46%) data are fit to Eq. 4. Better results are obtained iil/^o is allowed to be non-zero: 
^o(^ ~ 41%)/a = 50(20) and uriz ^ 41%) = 0.72(7). A value of ~ 0.7, together with 
zs=i/2 = Zp, is also obtained in a scaling analysis of our numerical data. We note that this value 
lies below ut = 0.90(5) for the S = 5/2 system Rb2(Mn,Mg)F4 [[T^. Fits to the data at the 
highest concentration give ut = 0.70(7) and ^o/a = 7(1). 

The percolation threshold Zp 40.7% is for the NN square lattice, and a non-zero frustrating 
next-NN (NNN) exchange could, in principle, shift Zp and lead to low-temperature spin-glass 
physics. The NNN exchange in La2Cu04 is expected to be frustrating and may be as large as 
8% J [ITSl], which would lead to a small reduction of the ordered moment, but should not shift the 



critical point [ [T0| ] . The dominant further-neighbor interaction may actually be a ring exchange 



among the four Cu sites on a square [POQ, which cannot extend connectivity beyond the NN 
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percolation threshold. In any case, further-neighbor interactions do not seem to noticeably 
affect the nature of the transition in La2Cui_2(Zn,Mg)204, , ,25=1/2 = Zp within the uncertainty 
of our experiment, and we find no evidence of spin-glass behavior. 

The correction terms to Eq. 1 can be expected to change as a result of the observed evolution 
of the low-temperature structure. Above z ~ 25%, Neel order occurs in a low-temperature 
tetragonal phase. Nevertheless, we find that Tn{z) evolves smoothly. Interestingly, our QMC 
data for Eq. 1 indicate that, as for undoped La2Cu04 WJ\, T^iz) corresponds very well to the 
temperature at which ^20 /O' = 100, as shown in Fig. 2B7This suggests that any changes in the 
correction terms in the full spin Hamiltonian must be very subtle. 

Assuming that quantum fluctuations do not alter the classical T = correlation length 
exponent, z/q = = 4/3 [|]], our result of ut ~ 0.7 suggests a crossover exponent $ = 
uq/ut ~ 1.9, which is larger than the values $ = 1.43 to 1.7 predicted theoretically |3T| ]. 



On the other hand, if z = Zp is a quantum critical point, it would suggest a dynamical critical 
exponent of 2; = I/ut ~ 1-4. Series expansion and numerical approaches for the S* = 1/2 
bond percolation problem, which is closely related to the site percolation problem of the present 
study, predict z = 1.7 to 2 [B^, whereas recent QMC work for the site-diluted 5=1/2 SLHAF 
arrives at 2; ^ 2.5 [P^. Our value is lower than these results, but also lies significantly above 
z = 1 for the Lorentz-mvariant QNLaM []2T1], implying that the QNLcrM fixed point is unstable 
to randomness. 

Recent theoretical work for the randomly diluted S = 1/2 SLHAF has led to some ex- 
act results in the dilute limit of quantum impurities and to interesting new predictions at 
higher concentrations [ |T2| , [T6| ]. It has been argued that the presence of impurities should lead 
to localized spin excitations and to the breakdown of the classical hydrodynamic description 
of excitations in terms of spin waves above a characteristic length scale [Jig]. In this picture, 
the spin-wave velocity c{z) is not a well-defined quantity, but there is no instability toward a 
disordered phase. Unlike dynamic observables, static properties such as the staggered magne- 
tization and T/v remain well defined throughout the ordered phase. Because the spin-stiffness 
also remains well defined in this theory, one continues to expect the correlation length to have 
the low-temperature form ^ ~ g^-^ps/T f^^. p^j-g system. We have been able to determine 
^(z,T) over a very wide range of impurity concentrations and temperature, and arrive at the 
unexpected conclusion that Eq. 2, and especially the heuristic crossover form Eq. 3, provide an 
excellent description of our data. 
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Table 1: Spin stiffness and spin- wave velocity extracted from fits of numerical data to Eq. 3. 
For z < 35%, ut was fixed to be 1. At 2; = 35%, a better fit was obtained using ut = 0.88(1). 
Equation 5 is modified by substitution of 1 + 2 for 1/Poo, as discussed in the text. Theory for 
c{z) is given by c(z)/c(0) = A{z)(l + z/2), as described in the text. The values for 27rps(0) 
and c(0) are from [^. 

z 2-Kps{z)/J Modified Eq. 5 c{z)/J/a Theory 

1.18(1) 1.13 1.33(3) 1.66 

0.08 0.79(1) 0.71 1.35(8) 1.31 

0.20 0.31(1) 0.28 1.43(9) 0.79 

0.31 0.063(1) 0.065 1.03(2) 0.34 

0.35 0.026(1) 0.025 0.69(1) 0.18 




Figure 1: Schematic of finite-sized sections of the infinite square lattice with random site 
dilution levels well below (31%), just below (40.7%), and above (45%) the percolation threshold 
Zp ~ 40.725% [|^]. Sites on the infinite cluster are shown in red, sites on finite disconnected 
clusters in blue, and diluents in white. The inset is a close-up view for z = 40.7%, showing the 
role magnetic Cu and non-magnetic Zn/Mg ions play in the experimental system. 
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Figure 2: (A) ( 1 ,0,0) magnetic peak intensity from neutron diffraction (orthorhombic notation) . 
A temperature-independent background has been subtracted. Lines represent fits for the mag- 
netic order parameter squared, ~ (T/v — T)^^ with /3 ^ 0.30, assuming a Gaussian distribution 
of Tat (typically ^ 4 K) to describe the rounding due to the small inhomogeneities present in 
the large samples used. (B) Tn from neutron diffraction (crystal) and magnetometry (crystal 
and powder sample). Above z ^ 20%, Tn{z) deviates from an extrapolated line, approaching 
zero at the percolation threshold Zp. Lines correspond to constant correlation lengths ^20/0- 
= 25 and 100 from Monte Carlo simulations of the randomly diluted S = 1/2 NN SLHAF; 
dashed regions are extrapolated from higher temperature. (C) Staggered moment per Cu atom, 
normalized by the value for the pure system. The neutron data are consistent with previous 
NQR results [|^]. The data are well described by a power law with zs=i/2 = Zp and exponent 
P^fj- = 0.45(3). Also shown is recent QMC and theory ^] for S" = 1/2, as well as the 
classical (S 00) result [[T5|]. 
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Figure 3: (A) Representative equal-time structure factor data. The solid red line represents 
a fit to the data, as discussed in the text, and the dashed black line indicates the instrumental 
resolution. This measurement was carried out in two-axis mode with 30.5 meV incident neu- 
tron energy and horizontal coUimations of 40'-27.5'-sample-23.7'. (B) Spin-spin correlation 
lengths in units of the lattice constant. Colored symbols represent results from neutron scat- 
tering measurements of La2Cui_2(Zn,Mg)204; black symbols represent Monte Carlo data for 
z = 8, 20, 31, 35, 41, and 46%. No adjustable parameters were used in the comparison. Exper- 
imental and numerical results for z = are from [IFtL ^ and the solid line is Eq. 2. Dashed 
lines are fits to Eqs. 3 and 4, as discussed in the text. 
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